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SUMMARY 

The linear aeroelastic equations for one curved blade of a vertical 
axis wind turbine are presented in state vector form. An attractive method, 
based on a simple integrating matrix scheme together with the transfer ma- 
trix idea, is proposed as a convenient way of solving the associated eigen- 
value problem for general support conditions. 


INTRODUCTION 

The troposkien shape has been frequently invoked in structural mod- 
eling of vertical axis wind turbine (VAWT) blades [1] . Although it is free 
of bending stress in the equilibrium position determined by a constant an- 
gular velocity, tests have indicated that under certain conditions serious 
vibrations about the original shape may occur. They Involve both chordwise 
and flatwise bending as well as torsion of the blade. 

Few theoretical works have been published on the determination of 
natural frequencies and vibration mode shapes of such special but neverthe- 
less important kinds of blade geometry [2-4]. The tendency among investiga- 
tors has been to deal with a more general geometry, which brings perhaps un- 
necessary complexity into the problem and requires the use of strictly nu- 
merical procedures for obtaining reliable results. Futhermore, the shape of 
actual VAWT blades has been very close to the troposkien; good approxima- 
tions like those suggested by Reis & Blackwell [5] are commonly used in prac- 
tice. One can recall that the complicated flow field generated by VAWT al- 
ready adds an important source of inaccuracy to the aerodynamic load esti- 
mation and consequently to problems such as flutter speed determination. 
Therefore it is important for preliminary design purposes, to develop semi- 
analytical methods based on the troposkien idealization, which will enable 
quick analyses and general understanding of the conditions that may affect 
blade stability. The intention of this paper is to present an approach with 
the aforementioned characteristics. 
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"OPTIMUM" LINEARIZED EQUATIONS OF MOTION 


The general second-degree nonlinear aeroelastic equations of motion 
for a slendei^ nonuniform and extensible blade were presented by Kaza & Kv- 
artenik [6], For physical understanding of the aeroelastic phenomena, how- 
ever, such degree of complexity does not seem necessary. It would be more 
interesting to seek a set of linear equations that allows the study of the 
relative importance of various parameters on the solution. These parameters 
include CG offset, shear deformations, rotary inertia and different kinds 
of blade support, along with the possibility of extending the model to a 
more realistic one containing the influence of tower and guy-cable stiff- 
ness as well as the other blades. 

When gravity effects are neglected, the troposkien can be seen as 
the plane curve described by a light rope rotating at constant angular veloc- 
ity. Its shape is known in closed analytic form [IJ. Working In the same 
fashion as Ashley [4], small perturbations from this equilibrium shape are 
assumed. After adding torsion to chordwise and flatwise bending and by as- 
suming simple harmonic motion in time, the final order of the governing sys- 
tem of differential equations is found to be 12. However, the choice of the 
dependent variables remains open. If, as in reference 4, the three cartesian 
components of the linear displacement related to the bending, the rotation 
related to the torsion and the incremental tension are selected, a set of 
coupled linear partial differential equations subject to one geometric con- 
straint equation is obtained. Asymptotic solutions may be useful, but only 
simplified versions of the problem could be carried out without a great deal 
of difficulty. An integrating matrix approach, like the one suggested by 
Hunter [7], could be employed in order to solve the eigenvalue problem for 
the free vibration analysis. The choice of the spatial independent variable 
as the distance along the blade introduces technical difficulties however, 
and the method fails. A more ’’natural” spatial independent variable is dis- 
covered to be the angle formed by the tangent to the undeformed shape and 
the axis of rotation. Even then, however, an unnecessary degree of algebra- 
ic difficulty is introduced. Without discarding the possibility of the so- 
lutions by the aforementioned means, the Reissner’s variational principle 
of the elasticity suggests that the best dependent variables are likely to 
be the 12x1 state vector of ’’natural” quantities connected with both equi- 
librium and deformation of the blade: bending moments, torque, shears, bend- 
ing slopes and displacements. The integrating matrix [7] together with the 
transfer matrix method [8] , is expected to constitute an efficient and pow- 
erful tool for dealing with the vibration of such beams subjected to gener- 
al boundary conditions. “ 
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AEROELASTIC ANALYSIS OF THE BLADE 


* 


Assuming small perturbations from the equilibrium shape is equiva- 
lent to constructing a linear field of deformations on a prestressed, un- 
twisted, plane troposkien-curved rod. The governing equations for a three- 
dimensionally curved rod under such a general initial stress configuration 
were obtained by different authors. Taking Nair & Hegemier*s derivation [9], 
based on the principle of virtual work, neglecting cross section warping 
and introducing the conditions defined by the troposkien hypothesis, one 
obtains a set of 12 linear first-order ordinary differential equations with 
nonconstant coefficients. These govern the 12 perturbation quantities cho- 
sen to be the dependent variables: 
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where Figure-1 together with the list below, define the nomenclature used: 

( )’= differentiation with respect to s; s= spatial independent variable, 
taken as distance along the locus of the section shear centers (elastic ax- 
is) ; Qj^_ 2 ~ internal shear components; T= tension along the blade; 
internal moment components; u,v,w= elastic axis displacements; ro- 

tations about the elastic axis (a= angle of attack); k = local Initial cur- 
vature, given by the troposkien; T = local initial tension, given by the 
troposkien; e= CG offset, constant^along the blade; El ,E1 = bending ri- 


gidities; EA= longitudinal rigidity; GJ= torsional rigflity; GA= effective 
shear rigidity; 2-3°^ external loads (forces and moments) of both 

inertial and aerodynamic origin. 

It became evident during the development that some extensibility 
should be allowed in this formulation in order to get the most rational set 
of equations. If inextensibility is assumed by simply letting EA become in- 
finitely large, important information is lost. Shear deformation is also in- 
cluded for completeness. 

The analysis proceeds by grouping the dependent variables in a 12x1 
state vector y, which contains generalized internal forces y^ and general- 
ized displacements y . The external loads are also collected in a 12x1 vec- 
tor p: _ 

y-hp|yDt - «i ‘>2 “ 2 ’■ ”3l“>‘i’" *2 ” “I 
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where f = f'' '^+f'’ respectively its aerodynamic and inertial components. 
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Hamilton's principle is involked to evaluate f as function of the 
generalized displacements. By taking both the translational and rotational 
parts of the kinetic energy, coriolis and rotary inertia effects are auto- 
matically included. After the rigid blade inertial load is subtracted from 
the equations, the following approximations appear to be significant to 
bring all inertial terms into a convenient form: 

« J'JpT\^d^dT\ C3 J n^)d?dn 


Here coordinates ^ and n are defined in Figure-1, and p is the blade mass . 
per unit of volume. 

As a first approximation, any reasonable aerodynamic theory is ade- 
quate to complete the aeroelastic problem formulation. The quasisteady strip 
theory can be employed without a great deal of difficulty. The author cites 
Kaza & Kvaternik^s derivation as appropriate [6]. 

Simple harmonic motion is assumed in order to eliminate the time de- 
pendence: this approach is known to be suitable for analyzing both free vi- 
bration and flutter. It follows that the inertial (J) and aerodynamic (A) 
loads can be written as: 


P 


[ J + A ] y = 



where and A^^^ are 6x6 submatrices dependent on the harmonic frequency. 
In Appendix A the nondimensional forms of the equations are presented for 
the free vibration problem in particular. 
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SOLUTION OF THE AEROELASTIC PROBLEM 


The main objective of this paper is to propose a convenient method 
for solving the foregoing aeroelastic problem. Numerical results are expect- 
ed to be published in a later work. Therefore, the author will now de- 
scribe how to obtain the eigenvalue problem associated with either free vi- 
bration or flutter analysis. 

Lehman [lO] presents a modified version of an integrating matrix 
scheme apparently first used by Hunter [7] which can be adapted for obtain- 
ing transfer matrices relating state vectors at different stations of a one- 
dimensional elastic structure. Appendix B gives a detailed derivation of the 
global trcinsfer matrix relating state vectors at the two ends of the 
blade: ° n 

^n = «o ^o 


. The enforcement of any physically possible set of boundary condi- 
tions at y and y will lead to the determination of the corresponding eigen- 
frequencies (as well as airspeed in the case of flutter). Such eigenf requen- 
cies are related to the complex roots of the so called frequency determinant 
associated with the particular set of selected boundary conditions. The ac- 
curacy of the method relies on the number of discretization points taken dur- 
ing the integration process, as well as on the basic polynomial from which 
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the integrating matrix is constructed (cf.ref.lO). The transfer niatrlx capa- 
bility for dealing with inhomogeneous boundary conditions, through the in- 
troduction of a new transfer matrix of very simple form expressing the lo-' 
'cal inhomogeneity, is what really makes this approach attractive. The con- 
struction of more elaborate models of VAOT, including for example guy wires 
and tower stiffness, seems possible without difficulty. Spring type bound- 
ary conditions are easily handled by a single matrix multiplication of the 
form: 
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,n+l 
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Here U is the new global transfer matrix and K is the spring transfer ma- 
trix. ° 


CONCLUDING REMARKS 

A semi-analytical method based on an integrating matrix scheme was 
proposed for obtaining the transfer matrix of a troposkien-curved blade. 
Such an approach is expected to be very efficient to study both free vibra- 
tions and flutter of VAOT. 
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APPENDIX A 
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NONDIMENSIONAL EQUATIONS FOR FREE VIBRATION IN STATE VECTOR FORM 

y’ = t s + J 1 y(i) 
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APPENDIX A 


where the new symbols are; a,= El /c^GA; a = [l-(x /r a„=EI /EI__ 

1 nn * a o j nn 

a^= EI^^/GJ; GA= shear rigidity; ( )= dimensionless quantities: (force) /V, 

(moment) /Vc, (displacement) /c; '( )= diagonal matrix; 1= identity matrix; 

T 

0= null matrix; ( ) = transpose matrix; (,)=zero element; x^=» dimension- 
less cross section static unbalance (in units of chord) ; r^*= dimensionless 

cross section radius of gyration (in units of chord) ; Ji= blade length; c= 
blade chord; m= blade total mass; V= vertical component of the initial ten- 
sion (constant for a gravity-free troposkien) ; f2= angular velocity of the 
VAOT; u= harmonic frequency; 4’~ angle defining the troposkien shape (Fig- 
ure-2); Y= shear coefficient; s= s/£; k = ; m= u/ft; -6 =sin x; C = cos x 

■i.= V^i”; sym= symmetric matrix. ° 
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APPENDIX B 


TRANSFER MATRIX DERIVATION WITH THE USE OF THE INTEGRATING MATRIX 
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In state vector form, the nondimensional aeroelastic equations can 
be written as: _ • 

y’=[S + J + Ajy = Zy ^ ^ ^ 
where S, J and A are 12x12 previously defined matrices and y is the 12x1 
dimensionless state vector of the dependent variables. When the integral 
operator L (cf.ref.ll) is applied, one gets: 


y = L Z y + k * (B-2) 

where k is a constant vector, dependent upon the boundary conditions .Pre- 
multiplying equation B-2 by a boundary condition matrix (cf .ref . 7) ,the 
state vector at one end is isolated since, by their intrinsic structure, 

D L = 0. Therefore y = k and B-2 can be re\7ritten as: 
o o _ 

y^ = ( I - L Z ) y (B-3) 

One notices that the term in parentheses is, by definition the inverse of 
the transfer matrix U. Hence, 


(B-4) 


y = U y = ( I - L Z y 
o *^o 

where I stands for the unit matrix. In particular, premultiplying U by an- 
other boundary condition matrix at the opposite end , yields: 

(B-5) 
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= D U y 
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n 


that is, the relation between state vectors at the two ends of the blade. 
The transfer matrix U can be seen as a product of intermediate transfer 
matrices relating conHitions at adjacent points along the integration 
path. The manipulation of most of these matrices is facilitated by their 
block structure and sparseness. 
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